Fast fluorescence dynamics in non-ratiometric calcium indicators 

M. Gersbach/ D. L. Boiko,^ * C. Niclass/ C. Petersen/ and E. Charboni 

Ecole Polytechnique Federale de Lausanne, Quantum Architecture Group, 1015, Lausanne, Switzerland 

^Centre Suisse d' Electronique et de Microtechnique SA, 2002, Neuchatel, Switzerland 

* Corresponding author: dmitri.hoiko® csem.ch 



00 

o 

O 
(N 

> 

O 



Oh 
I 

O 



c/3 
O 

Oh 



(N 
> 

in 

(N 

o 

00 

O 
> 

X 

?-( 



Compiled November 13, 2008 

The fluorescence decay of high-affinity non-ratiometric Ca^+ indicator Oregon Green BAPTA-1 (OGB-1) is analyzed 
with unprecedented temporal resolution in two-photon excitation regime. A triple exponential decay is shown to best fit 
the fluorescence dynamics of OGB-1. We provide a model for accurate measurements of the free Ca^+ concentration and 
dissociation constants of non-ratiometric calcium indicators. © 2008 Optical Society of America 
OCIS codes: 170.6280, 160.2540, 260.2510, 170.2520, 300.6280, 180.2520. 



Fluorescence Lifetime Imaging Microscopy (FLIM) is 
used to locally probe the chemical environment of fluo- 
rophores, e.g., ion concentration, pH, or oxygen content 
[1, 2]. To acquire time-resolved fluorescence images, the 
technique of Time-Correlated Single Photon Counting (TC- 
SPC), in combination with detectors exhibiting single-photon 
sensitivity [3], is commonly used. This technique enables 
the measurement of photon time-of-arrival distributions with 
very high accuracies, independently of instabilities in the 
excitation beam intensity. So far, temporal resolutions of a 
few hundred picoseconds were considered sufficient in bio- 
medical FLIM applications. 

In this letter, we analyze the fluorescence decay of the 
high-affinity Ca2+ indicator Oregon Green BAPTA-1 (OGB- 
1) under two-photon excitation conditions, using a TCSPC 
system of measured resolution 79ps based on CMOS Sin- 
gle Photon Avalanche Diode (SPAD) detector technology. 
Our measurements reveal a triple-exponential decay of OGB- 
1 fluorescence, which we show enables accurate measure- 
ments of Ca^+ concentration and dissociation constants of 
non-ratiometric fluorescent probes. We provide a comparison 
with previously reported data [1,4,5], which were acquired at 
resolution of 200ps. 

The SPAD-based TCSPC system is depicted in Fig.L Our 
SPADs [6] are integrated in a 32x32 array and incorporate 
on-chip high bandwidth I/O circuitry. The active region of a 
SPAD pixel [Fig. 1(b)] consists of a p+-n junction operating in 
the Geiger mode. Due to a small diameter (7 /im), our SPADs 
show an extremely low dark count rate (< 10 Hz at room tem- 
perature). The photon detection probability is 25 % at 500 
nm wavelength, while the dead time is 25 ns with negligible 
afterpulsing (< 0.1 %) [6]. 

Fluorescent molecules are excited in the two-photon 
regime [7], using a mode-locked Ti:Sapphire laser (MaiTai, 
Spectra Physics) emitting 100 fs optical pulses at 800 nm 
wavelength [Fig. 1(a)]. The attenuated beam with an average 
power of 9 mW is focused on a sample using a 20x micro- 
scope objective (XLUMPlanFL, Olympus), which also serves 
to collect the fluorescent emission. Fluorescence from OGB-1 
samples treated here (central emission wavelength 520 nm) is 
directed towards the SPAD by a dichroic beam splitter (DBS) 
and a filter (E650SP, Chroma Technology) for suppression 



of backscattered excitation pulses. Another 20x objective im- 
ages the emission spot onto the SPAD. The low quantum yield 
of the two-photon excitation fluorescence results in a count 
rate at the detector of 10 kHz for 80 MHz repetition rate 
of excitation pulses. Time discrimination is performed in a 
"reversed start-stop" configuration by measuring the time in- 
tervals from a photo-detection event at the SPAD to emission 
of a successive laser pulse. For precise timing of photon ar- 
rivals, we use a 6 GHz bandwidth oscilloscope (WaveMas- 
ter, LeCroy) incorporating a time-to-digital converter and en- 
abling computation of histograms of photon arrivals. 

Fig. 1(c) shows the Instrument Response Function (IRF) of 
the entire system at 400 nm wavelength (blue curve). It is 
recorded by measuring the hyper-Rayleigh scattering of 800 
nm wavelength pulses in a solution of colloidal gold parti- 
cles (G1652, Sigma- Aldrich) [8]. At the incident power of 90 
mW, the average count rate of the detector is just 600 Hz. The 
measured photon arrival time jitter of 79 ps (FWHM) is dom- 
inated by the SPADs time-response characteristics. The IRF 
is only slightly asymmetrical and assumes a Gaussian-curve 
approximation IRF— cxp(— i^/2o'j^]^p) (inset, cyan curve). 

Fluorescent samples were composed of 2 /il OGB-1 dye 
(O6806, Molecular probes) and 20 ^1 Ca:EGTA buffer solu- 
tions from the calibration kit (C3008MP, Molecular Probes) 
with quoted free Ca^+ concentrations in the range of 17 nM 
to 39 /iM. Several fluorescence decay curves measured in the 
two-photon excitation regime are shown in Fig. 1(c). The fluo- 
rescence lifetimes were obtained from the numerical analysis 
of these data. 

Thorough numerical fit must take into account the IRF 
of the system [5]. As opposed to conventional systems with 
strongly asymmetric IRF, the Gaussian-like IRF of our sys- 
tem assumes a simple analytical expression for the measured 
fluorescence decay. For a train of excitation pulses of period 
T, each term of a multi-exponential decay process reads 
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where a triple-exponential decay is assumed and index k in- 
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Fig. 1. (Color online) (a) Schematic of the experimental 
setup for fluorescence lifetime measurements, (b) Scanning 
electron microscope image of a SPAD. (c) Measured fluo- 
rescence decay of OGB-1 at various calcium concentrations 
(gray scale curves) and hyper-Rayleigh scattering from col- 
loidal gold particles (blue curve). The red curve highlights 
measured fluorescence decay in a 38 nM free Ca^+ concen- 
tration buffer and the numerical fits of the double-exponential 
(olive curve, x^=1.80) and triple-exponential (green curve, 
Xr=l-09) decay models. The inset shows a close-up of the 
initial interval of 600ps width and (d) shows the residues for 
the two models. 



dicates fast (/), intermediate (i) and slow (s) temporal com- 
ponents, Tfe is the fluorescence emission hfetime (Tf<Ti<Ts) 
and erf(z) = -^ L exp(— ^^)(i^ is the error function. Eq.(l) 
takes the periodic train of excitation pulses and response time 
jitter of the SPAD into account, such that our data in Fig.l do 
not require deconvolution processing. 

Analysis of the OGB-1 fluorescence reveals the best agree- 
ment with a triple exponential decay approximation. The data 
are modeled using the function AfIf+AiIi+{l~Af—Ai)Is 
with the fast If{t), intermediate /^(t) and slow Is{t) de- 
caying components (1) of normalized partial intensities Af, 
Ai and As=l~Af~Ai, respectively. Fig. 1(c) details a com- 
parison between the double-exponential (olive curve) and 
triple-exponential (green curve) models applied to fluores- 
cence from a 38 nM Ca^+ concentration sample (red points 
and curve). In both cases, the residues [Fig. 1(d)] do not show 
any bias that might be caused by the Gaussian curve approx- 
imation of the IRF, while the quality of the numerical fit 
is improved in case of the triple-exponential model, as also 
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Fig. 2. (Color online) Partial intensities ^4^, Ai and Af of the 
slow (rs=4. 18±0.01ns), intermediate (Ti=768±16ps) and fast 
(r/=188±6ps) decay components (a) and the mean lifetime 
T-iw&—J2k^kTk (b) of OGB-1 fluorescence as a function of the 
free Ca^+ concentration in comparison with the data from [4] 
and [5], Xr=0.96. The numerical fit to the model (2) (curves) 
yields i^D=195±9 nM, KDn=4±6 ^iM and n=9±3. (c): The 
ratio Aj[Ca^^]/y4s (points) and its variance-weighted mean 
yielding /•ir£)=180±llnM (dashed line). The inset shows the 
structure of OGB-1 molecule. 



confirmed by the reduced x^-values. Wilms et al [5] have 
made the same observation for OGB-1 fluorescence in the 
absence of Ca^+. However, it was attributed to contaminat- 
ing dye derivatives and a double-exponential model has been 
utilized for samples containing Ca'^^. Previously reported 
data [1,4,5] for OGB-1 fluorescence lifetimes thus assume 
double-exponential decay. We argue that at high Ca^+ con- 
centration, large amplitude of long-leaving component make 
it difficult to detect short lifetime components with small am- 
plitudes and the FLIM systems used in [1,4,5] exhibit insuf- 
ficient temporal resolution to enable observation of the triple- 
exponential behavior 

Using triple-exponential fluorescence decay (1) with the 
lifetimes independent of the calcium buffer [1,5], we ap- 
plied a global numerical analysis to our data, yielding the 
partial intensities Ak in function of free Ca^+ concentra- 
tion and the lifetimes r^ (Fig. 2). Temporal resolution of our 
system allows the short lifetime component (Ty~188±6ps) 
to be unambiguously resolved in the background of the 
intermediate- and long-living fluorescence (Ti~768±16ps 
and T<;~4.18±0.01ns). The partial intensities Af and Ai de- 
crease while the slow-component intensity As=l—Af—Ai 
increases with concentration [Ca^+]. 



We attribute, as usual, the short- and long-Ufetime com- 
ponents to, respectively, unbound and Ca^+-bound OGB-1 
with a simple 1:1 complex stoichiometry. As such and be- 
cause the intensities are defined by complex concentrations, 
we obtain that A/cx[D] and As(x[CnD]^Kjj^[Cn'^+][D] with 
[D] and [CaD] being the unbound and Ca-bound OGB-1 con- 
centrations and Kd the dissociation constant of the 1 : 1 com- 
plex. The ratio A/[Ca ^]/As thus yields an estimate of the 
dissociation constant. In Fig. 2 (c), a weighted average ratio 
180±llnM is in good agreement with the quoted (by manu- 
facturer) Kd of 170 nM. 

In Fig. 2, the long lifetime component As is in good agree- 
ment with the measurements of Agronskaia et al [4] and 
Wilms et al [5], but in [5], As<l in Ca^+ -saturated buffer 
and Kd is too high (^300 nM). These discrepancies are at- 
tributed to the dye impurities and to the difference in two- 
photon absorption cross-sections of Ca-bound and unbound 
OGB-1, which are said difficult to be quantified [5]. 

For the Ca-bound OGB-1, Lakowicz [1] reported t^ of 4 ns, 
which agrees well with our data. In [4], Tg varies in the range 
2.6-3.7 ns, while in [5], Ts=3.63 ns. For the Ca-free OGB-1, 
a single lifetime component Tf has been measured of 700 [1], 
290-420 [4] and 346 ps [5], respectively. These results cor- 
respond to the combined effect of the decay processes with 
the lifetimes Tj^l90ps and ri^770ps in our measurements, 
which may not have been resolved in previous studies. The 
experimental setup in [5] relies on a commercially available 
Photomultiplying Tube (PMT) with a quoted timing jitter of 
200 ps. Other data [1,4] are reported without timing resolu- 
tion of experimental setups. 

Due to the Ca-binding features of the octadentate chelator 
BAPTA [9], the 1:1 stoichiometry is usually attributed to the 
Ca-bound OGB-1 as well. However, the molecule of OGB-1 
has an asymmetric structure with respect to its BAPTA moi- 
ety (Fig. 2, inset). The asymmetric arrangement of the car- 
boxyl functional groups, benzol rings, fluorine and nitrogen 
atoms as well as the variations of electronic density from low 
(at hydrogen in carboxyl groups) to high (at benzol rings, F 
and O atoms) assume a polarity of the molecule and, as a 
consequence, weak dipole-dipole intermolecular forces. As 
a result of such interaction forces, several OGB-1 molecules 
might be coordinated to a Ca-bound OGB-1 forming thus a 
polymolecular association with a calcium:indicator molar ra- 
tio l:n. The intermediate lifetime component t^ in Fig. 2 then 
might be attributed to such polymolecular association, yield- 
ing A,oc[Cai/„D]=X^y"[D][Ca2+]i/n with Kon being the 
corresponding dissociation constant. Note that formation of a 
polymolecular structure is a very complex processes and re- 
quires thorough investigations but a non 1 : 1 stoichiometry of 
Ca:OGB-l has been reported in [10] and here it allows us to 
build an accurate model: 
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In the limit Kun^oo (double exponential decay), it agrees 
with the Hill equation, as opposed to the model in [5]. The 
numerical fit of data in Fig. 2 reports small binding affinity 
(KDn^4:fiM) and n=9 indicating the most probable form of 



polymolecular association with eight Ca-free molecules co- 
ordinated to the Ca-bound OGB-1. The dissociation constant 
K]j reported by the fit is 195 nM, in agreement with the av- 
erage of Af [Ca^^] /As (bottom panel). 

The triple exponential fluorescence decay in non- 
ratiometric Ca probes as a result of dye contaminations was 
suggested by Lakowicz et al [11] for the Calcium Green 
(CG-1). The CG-1 molecule has the same structure as OGB- 
1 but the two F atoms are replaced by CI, yielding [11] 
T/=50ps, Ti=450ps and rs=3.7ns. Interestingly, the signa- 
ture of non- 1 : 1 stoichiometry has been noticed for Ca-bound 
CG as well [12]. We find that at low Ca^+ concentration, 
the ratio Af[Ca^+]/As in [11] is a constant of 170±20 nM, 
yielding Kjj close to the quoted value of 190 nM. At high Ca 
concentration, the band pass of the system (2GHz modulation 
using a frequency-domain FLIM technique) was insufficient 
to accurately measure the small component Af. The data in 
[11] can thus be accurately interpreted in the framework of 
our model (2), without appealing to dye impurities. (K]j ob- 
tained in [11] using a conventional model [13] is 128nM). 

In summary, we have shown that high-temporal resolution 
measurements of the triple-exponential fluorescence decay 
of non-ratiometric Ca^+ indicators allow the free Ca^+ con- 
centration and dye dissociation constant to be measured pre- 
cisely. 
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